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We derive the dynamical mean-field theory equations for transport in an ordered charge-density- 
wave phase on a bipartite lattice. The formalism is applied to the spinless Falicov-Kimball model on 
a hypercubic lattice at half filling. We determine the many-body density of states, the dc charge and 
heat conductivities, and the optical conductivity. Vertex corrections continue to vanish within the 
ordered phase, but the density of states and the transport coefficients show anomalous behavior due 
to the rapid development of thermally activated subgap states. We also examine the optical sum 
rule and sum rules for the first three moments of the Green's functions within the ordered phase 
and see that the total optical spectral weight in the ordered phase either decreases or increases 
depending on the strength of the interactions. 
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I. INTRODUCTION 

Dynamical mean-field theory was introduced almost 
two decades ago by Brandt and Mielsch 1 , who solved 
for the transition temperature into a charge-density-wave 
(CDW) phase of the spinless Falicov-Kimball model at 
half filling. This work appeared shortly after the idea 
of examining strongly correlated electrons in the limit 
of infinite dimensions was introduced^. Since then, the 
field of DMFT has emerged as one of the most pow- 
erful nonperturbative techniques for solving the many- 
body problem. While results for many properties exist 
in the homogeneous (unordered) phased, there has been 
little work in examining the properties of the ordered 
phase. Brandt and Mielsch worked out the formalism for 
calculating ordered-phase Green's functions^, the order 
parameter was shown to display anomalous behavior at 
weak coupling 5 ' 6 , and higher-period ordered phases have 
been examined on the Bethe lattice 7 . But, surprisingly, 
there has been no work on the transport properties in 
the ordered phase. Indeed, it is interesting to compare 
how transport varies in the homogeneous phase versus 
the ordered phase. At weak coupling, we anticipate the 
gap formation of the CDW to greatly suppress the dc 
transport, while at strong coupling it may be a much 
milder correction to the Mott-insulating behavior. What 
is more interesting is to examine the temperature depen- 
dence. For example, in systems that are metallic at high 
temperature, the many-body DOS in the CDW phase 
develops strong temperature dependence (with increas- 
ing T) as the CDW gap region fills in due to thermal 
excitations, until gap closure is complete at the transi- 
tion temperature. But unlike the well-known supercon- 
ducting case, where subgap states tend not to form and 
the gap is simply reduced in size as T increases, here we 
have a rapid development of subgap states, even though 
the CDW order parameter remains nonzero. These sub- 
gap states should produce anomalous behavior in the 



low-T transport, and indeed we find this is so but the 
quantitative behavior is not that different from expo- 
nential activation of the transport. We anticipate our 
results should be relevant to different experimental sys- 
tems that display charge-density-wave order, especially 
in compounds which are three-dimensional like 8 BaBi03 
and Bai-zK^BiOs. 

This contribution is organized as follows: In Section 
II, we present the formalism for DMFT in the ordered 
phase including the techniques needed to determine the 
optical conductivity and the dc transport. We also de- 
termine moment sum rules for the Green's functions in 
the ordered phase. In Section III, we apply the formal- 
ism to numerical solutions of the Falicov-Kimball model 
at half filling and show how the transport behaves in the 
ordered phases. Conclusions and a discussion follow in 
Section IV. 



II. FORMALISM FOR THE ORDERED PHASE 

The Falicov-Kimball model 9 was introduced in 1969 as 
a model for metal-insulator transitions in rare-earth com- 
pounds and transition-metal oxides. The spinless version 
is arguably the simplest many-body problem that never- 
theless possesses rich physics including the Mott transi- 
tion, order-disorder phase transitions, and phase sepa- 
ration (for a review see Ref. M). It involves two kinds 
of electrons: mobile conduction electrons whose creation 
and destruction operators are d\ and d t at site i; and 
localized electrons whose creation and destruction oper- 
ators are // and f t at site i. The Falicov-Kimball Hamil- 
tonian can be represented in terms of a local operator 
and a hopping operator as follows 

i ij 
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where Uj is the hopping matrix and 

Hi = Uh ld n l} - HdTkd ~ M/» 



(2) 



is the local Hamiltonian with the number operators given 
by ha, = d\di and h i} = fj f v 

If the lattice can be divided into two sublattices, and 
the hopping is nonzero only between the two sublattices 
(i. e., there is no hopping within either sublattice), then 
the lattice is called a bipartite lattice, and it has nest- 
ing at half filling in the noninteracting system, which 
implies the fermi surface in the Brillouin zone has flat re- 
gions that are connected by the zone-diagonal wavevec- 
tor Q — (tt,tv,.. .). Nesting promotes the formation of 
a CDW with the average filling of the electrons being 
uniform on each sublattice, but changing from one sub- 
lattice to another. This is often called the checkerboard 
or chessboard CDW, and is the ordered phase that we 
will examine in detail in this work. 
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FIG. 1: Schematic illustration of the transition from a homo- 
geneous phase to the bipartite CDW phase. The hopping is 
between nearest neighbors, which corresponds to the neigh- 
boring points in the horizontal and vertical directions. 

In order to develop the formalism to determine the 
Green's functions and transport in the ordered CDW 
phase, we need to introduce some notation that will help 
clarify how the ordered phase is determined. It is con- 
venient to supplement the lattice site index, which we 
had been calling i, by a double index (i, a), where i 
runs over all of the lattice sites of one of the sublattices, 
and the label a = A or B denotes the sublattice (see 
Fig. [TJ we are assuming for simplicity that the two sub- 
lattices have an equal number of lattice sites as they do 
on the infinite-dimensional hypercubic lattice or on the 
infinite-coordination- number Bethe lattice). We rewrite 
the Hamiltonian from Eq. ((TJ) as 



ia ijab 



(3) 



with the local Hamiltonian satisfying 

^ = « r /iK|-^/i (4) 
in this notation, the bipartite lattice condition is sim- 



j.AA _ J.BB 

chemical potentials for the two different sublattices at 



ply that tfj = tfj = 0. We have introduced different 



the moment. This is convenient for computations, be- 
cause it allows us to work with a fixed order parameter, 
rather than iterating the DMFT equations to determine 
the order parameter (which is subject to critical slowing 
down near T c ). Of course, the equilibrium solution oc- 
curs when the chemical potential is uniform throughout 
the system (fj,^ = /zf and fj,4 = fx?). 

Our starting point is to find the set of equations satis- 
fied by the lattice Green's function. The Green's function 
is defined to be 



GtKr) 



Tr 



T T e-^d la (r)d] b (0)} /Z, (5) 



where r is the imaginary time, the time dependence of the 
destruction operator is written in the Heisenberg repre- 
sentation {d ia (r) — exp[r7i]d ia exp[— rti]}, and Z is the 
partition function Z — Trexp[— /3n\, with (3 — 1/T the 
inverse temperature. The symbol T T is the time-ordering 
operator, which orders the times so that earlier times 
appear to the right. 

One way to calculate the Green's function is to use 
an equation of motion technique 11 , where the derivative 
with respect to imaginary time is taken and a differential 
equation is found for the Green's function. In DMFT, 
this procedure is carried out for the impurity problem in 
a time-dependent field, and the field is adjusted so that 
the impurity Green's function is equal to the local lat- 
tice Green's function. In addition, we need to define the 
self-energy via Dyson's equation in order to complete the 
iterative DMFT loop needed to solve the full problem. 
Finally, an analytic continuation from the imaginary axis 
to the real axis is performed to calculate dynamical prop- 
erties. These techniques are all well known and have been 
established in the lit er at ur e^ i 1 ' 1 1 , so we provide just a 
schematic approach to the derivation, highlighting some 
key formulas along the way. 

The Dyson equation, which can be thought of as defin- 
ing the self-energy is 



E 



[( W + MSacSu - esth + *S c ]G?jV) = M«*> ( 6 ) 



with uj the real frequency. In the case of nearest-neighbor 
hopping on an infinite-dimensional hypercubic lattice, we 
have that the band structure satisfies = — ^ . exp[ik • 



(Ri 



D 



liniD^oo t* cosfca/v^D, where 



we scaled 2 the nearest neighbor hopping matrix element 
by t = t* /2>/D (we will use t* — 1 as our energy unit). 
In addition, the self-energy is locate 2 - 



(7) 



which further simplifies the Dyson equation. It is sim- 
pler to transform from real space to momentum space to 
solve the Dyson equation. But we do not assume that 
the Green's function is completely translation invariant, 
instead, we assume only that there is translation invari- 
ance within each of the sublattices. Then the momentum 
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representation of the Dyson equation in Eq. Q with the 
local self-energy in Eq. ([7]) becomes 



Gfc(w) = [z(uj) -tk] 1 , 



(8) 



where z(lu) and the hopping term are represented by 2 x 2 
matrices 
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(9) 



Substituting Eq. © into Eq. ([8|) and taking the matrix 
inverse yields the following formulas for the momentum- 
dependent Green's functions on the lattice 



G£ B (u) = G%*(u) 



BA i 



(10) 

(11) 
(12) 



Z(u>) = ^[u; + (4- Z a (uj)}[u> + /zf - E*( W )], (13) 

which agree with those of Brandt and Mielsch^ even 
though our notation is somewhat different from theirs. 
The local Green's functions on each sublattice then sat- 
isfy 



with Z defined by 



G aa (io) = 



Fm(u), (14) 



where the a sublattice is different from the b sublattice 
and F z (cu) is the Hilbert transform 



(15) 



F z (u) = / dep(e) 



Z(tu) 



The function p(e) is the noninteracting density of states, 
which is p(e) = exp(-e 2 /t* 2 )/t*y/ir for the infinite- 
dimensional hypercubic lattice (as discussed above, we 
take t* = 1). 

In the DMFT solution, we need to map the lattice 
problem onto a local (impurity) problem in a time- 
dependent field that is adjusted to make the impurity 
Green's function equal to the local Green's function of 
the lattice. Here, we have two different local Green's 
functions, one on the A sublattice and one on the B 
sublattice; hence we will need two time dependent fields 
and two impurity problems to solve in order to complete 
the DMFT mapping. We call the dynamical mean fields 
\ a (uj) for each sublattice. Then the solution of the im- 
purity problem is straightforward and is summarized by 



the following set of equations 

Gg(w) = [(? 00 (w)- 1 +I! B (w)]~ 1 
1 



G aa {uo) 



(1 - n)) 
u + pf d - X a (u) 



£ a (w) = lo + p a d - A<» - G aa (Lo)-\ 



(16) 
(17) 



(18) 
(19) 



where we must solve these equations for each of the sub- 
lattices a = A and a = B. 

The DMFT algorithm for a fixed value of the order pa- 



A ' n? 



rameter starts by choosing and rij such that n* -r 
is fixed to the total /-electron filling (the order parameter 



is Anf = 



7 



n?), and choosing p A = fi d . With those 



fixed quantities, we propose a guess for the self-energy on 
each sublattice, and then compute the local Green's func- 
tion on the real axis from Eqs. lfl3|) and (fT4|) . Then we 
extract the dynamical mean field on each sublattice from 
Eqs. (TT5]) and (fT7|) . then find the local Green's function 
for the impurity from Eq. (|18|) and the new self-energy 
from Eq. (fT9| . This loop is repeated until the Green's 
functions converge. Then one can calculate the filling 
of the d-electrons and adjust them until they match the 
target filling. But this procedure is not yet complete, be- 
cause we need to determine the correct equilibrium order 
parameter n A — at the given temperature. To find 
this, it is actually more convenient to perform the calcu- 
lations precisely as described above, but on the imaginary 
frequency axis, where u is replaced by iuj n = iirT(2n + 1) 
the fermionic Matsubara frequencies. Then we calculate 
the chemical potential for the /-electrons on each sublat- 
tice via 



U 
2" 



Tin- 



1 
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T^ln[l (20) 



and adjust the order parameter until the two chemical 
potentials are equal, which is required for the equilibrium 
solution. Then, when we calculate the Green's functions 
on the real axis, the chemical potentials and fillings are 
all already known, so they do not need to be adjusted 
during the calculation. 

This algorithm is much more efficient than an algo- 
rithm that starts with a fixed chemical potential for the 
/-electrons and iterates to determine n°j on the imagi- 
nary frequency axis. This is because the latter suffers 
from critical slowing down, and becomes quite inefficient 
near the critical temperature, whereas the calculations 
with the fixed order parameter converge quite rapidly re- 
gardless of how close one is to the critical point. 

When the DOS is calculated for each sublattice in 
the ordered phase, one finds interesting temperature de- 
pendence of the subgap states as a function of T. It 
is illustrative to discuss these evolutions in terms of 



4 



the moments of the local interacting DOS. It is well 
known, that in the homogeneous phase, the integral of 
A(uS) — —JmG(u))/n is equal to 1. But there are also 
exact results known for higher moments as wel l 13 i 14 . In 
particular, because the moments are derived from oper- 
ator identities, they continue to hold whether in the or- 
dered phase or not. So we learn the following identities 
immediately: 
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dujA a (uj) 
dujujA a (uj) 

dujL0 2 A a (Lu) 



(21) 



(22) 



+ i xf-2U i x a d n a f + U'n a f . (23) 



We have checked these moments versus our numerical 
calculations of the Green's functions on the real axis and 
they all agree to high accuracy for all temperatures that 
we consider. Note that at half filling, we have p? d = U/2, 
so the first moment vanishes in the homogeneous phase. 
As the system orders, the first moment on one sublatticc 
becomes negative, and the first moment on the other sub- 
lattice becomes positive, which indicates that the quan- 
tum states are shifting in response to the ordering. In 
particular, this redistribution of states causes the average 
kinetic energy to evolve more strongly with temperature 
in the ordered phase, but its evolution is anomalous, and 
cannot be predicted by any simple reasoning about how 
the states evolve (see below). The evolution of the av- 
erage kinetic energy plays an important role in the total 
spectral weight for the optical conductivity. 

At T = 0, the order parameter goes to 1, so there is one 
sublattice (let us say the A sublattice) which has all the 
/-electrons. Hence — 1 and = 0. In this case, the 
analysis for the Green's function simplifies. In particular, 
only one term in Eq. (|18[) survives on each sublattice and 
we immediately find T, A — U and E s = 0. Plugging 
these results into the remaining formulas for the DMFT 
algorithm then yields an analytic formula for the ordered 
phase DOS 



Aa,b(u) 



1 



ImG 



AA,BB 



Re 





(24) 



where the top sign is for the A sublattice (with a diver- 
gence of the DOS at u = U /2) and the bottom sign is for 
the B sublattice (with a divergence at ui = —U/2); the 
formula is restricted to half filling where fi d = fi d = U/2. 
Note that the two DOS on each sublattice are mirror im- 
ages of each other and that each sublattice has weight 
for positive and negative frequency, but the band that 
does not have the singularity (lower band for sublattice 
A and upper band for sublattice B) has shrinking spectral 
weight as U becomes large, because the mobile electrons 
avoid the sites with the localized electrons for large U . 
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FIG. 2: (Color online) DOS at T = for the CDW-ordered 
phase on a hypercubic lattice. Panel (a) is for the A sublattice 
and panel (b) is for the B sublattice. Four cases are plotted: 
U = 0.5 which is a strongly correlated metal; U = 1, where 
a dip develops in the normal-state DOS at the chemical po- 
tential; U = 1.5, which is a near-critical Mott insulator; and 
U = 2.5, which is a moderate-size-gap Mott insulator. The 
T = gap in the DOS is always equal to U in the ordered 
CDW phase. 



Note further, that unlike the Mott insulator, where the 
DOS vanishes only at the chemical potential on a hyper- 
cubic lattice, a real gap develops here of magnitude U 
at T = 0. In Fig. [21 we show the DOS at zero temper- 
ature for four values of U. Panel (a) plots the DOS on 
the A sublattice and panel (b) plots the DOS on the B 
sublattice. One can see that the shape of the DOS is 
qualitatively similar for all cases, but the size of the gap 
grows with U . 

What is more interesting is to examine the tempera- 
ture evolution of the DOS in these different cases. Indeed, 
the system develops substantial subgap DOS that is ther- 
mally excited within the ordered phase (the order param- 
eter is determined by the difference in localized electron 
filling on the two sublattices). In Fig. [3] (a), we plot the 
DOS for the strongly correlated metal at U = 0.5. The 
fill in of the subgap states is quite rapid with T as we 
increase up to T c — 0.0336. Similar behavior is also ob- 
served for U = 1 with T c — 0.0615 which has a dip in the 
DOS in the normal state [Fig. [3] (b)]. 



5 



10 



< 

co 
O 
Q 



■ (a) 


— T=0.01 - 


— T=0.02 




T=0.03 


_ 


— T=0.04 

I. 


^2 





1 2 
Frequency [t*] 



10 



< 

CO 

O 
Q 



" (b) 






T=0.02 - 








T=0.04 








-T=0.06 








T=0.07 


















-2 -1 





1 


2 3 i 



Frequency [t*] 



10 



< 

CO 

O 
Q 



" (a) 




T=0.03- 






T=0.05 






-T=0.07. 






T=0.08 









1 2 
Frequency [t*] 



< 

CO 

O 
Q 



1 2 
Frequency [t*] 




FIG. 3: (Color online) DOS on the A sublattice for various 
T values in the CDW-ordered phase on a hypercubic lattice 
with (a) U = 0.5 and (b) (7 = 1. The DOS on the B sublattice 
is a mirror reflection of these results about the plane uj — 0. 



FIG. 4: (Color online) DOS on the A sublattice for various 
T values in the CDW-ordered phase on a hypercubic lattice 
with (a) U = 1.5 and (b) U = 2.5. 



The Mott insulating phases also illustrate interesting 
behavior. In particular, the subgap states develop pri- 
marily within the upper and lower Hubbard bands (al- 
though on the hypercubic lattice, the Mott insulator has 
only a pseudogap with the DOS strictly vanishing only at 
co = 0). We illustrate this behavior in Figs. [4] (a) an d (b). 
The transition temperatures are T c = 0.0747 for U = 1.5 
and T c = 0.0724 for U = 2.5. Note how the subgap DOS 
develop closer to the Mott band edge than they do to 
the CDW band edge, which implies they should have an 
effect on the transport at low T. 

In all cases, the DOS satisfy the three sum rules for the 
first three moments to essentially machine accuracy — our 
actual accuracy is determined by the step size we use for 
the real frequency axis in calculating the DOS and then 
integrating it over all frequency to obtain the numerical 
moments. 

Now we develop the formalism for transport in the 
CDW phase. The linear response optical conductivity is 
determined (via the Kubo-Greenwood formul a 15 i 16 ) by 
the imaginary part of the analytic continuation of the 
current-current correlation function to the real axis, 

= -ImUjj(uj), (25) 



with the number current operators defined by 

j = i J2 $ {R* - Rjb)$Jib , (26) 

ijab 
f) f ab 

^=E^4(fe)4(fe). (27) 

abk 

The procedure to determine the current-current corre- 
lation function is a standard one so we only sketch the 
derivation briefly. We start from the imaginary time for- 
mula for the current-current correlation function 

n jj (r-r') = {T r j(r)j(T')), (28) 

where the angle brackets denote a trace over all states 
weighted by the statistical operator (density matrix) at 
the given temperature and the current operators are rep- 
resented in the Heisenberg representation with respect 
to the equilibrium Hamiltonian (because this is a linear- 
response calculation). We then perform a Fourier trans- 
formation to go from imaginary time to Matsubara fre- 
quencies, and then perform an analytic continuation from 
the imaginary frequency axis to the real frequency axis. 

The Fourier transform of the current-current correla- 
tion function defined in Eq. (|28p can be represented as a 
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FIG. 5: Bethe-Salpeter equation for the generalized polariza- 
tion. 



summation over Matsubara frequencies 



(29) 



where we introduced the shorthand notation Il m m+ i = 
TL(iu) m ,iui m + ivi) for the dependence on the fermionic 
iuj m — iirT(2m + 1) and bosonic ivi — ilixTl Matsub- 
ara frequencies (m and / are integers). In the CDW 
phase, the graphic depiction of the Bethe-Salpeter equa- 
tion for the generalized polarization n mjm +„ is plotted in 
Fig. O where the solid oval depicts the current operator 
using the same sublattice indices as we have used before 
(the current operator connects the two sublattices), the 
solid lines are Green's functions, and the cross hatched 
object is the total (reducible) charge vertex. The cur- 
rent operator vertex contains the factor dt\~jdk a which 
is an odd function of the wavevector. Since the band 
structure and the Green's functions are even functions 
of the wavevector, any summation over momentum that 
contains one current vertex and any number of Green's 
functions will vanish. Now, in infinite dimensions, the ir- 
reducible charge vertex (which enters the Bethe-Salpeter 
equation for the total charge vertex) is local and hence 
momentum independent, so the second term in Fig. 
vanishes, just like it did in the homogeneous phased. 
We thereby conclude that the optical conductivity is con- 
structed only by the bare bubble in Fig.O Then, the full 



k,ia<„ 



^-m,m+l — /, j 
abed 




k, ico„+iv. 



k, icq,,, k, iC0» 

\ ,•• 1 A, 



=J 




J + J 



B 




J 



+J 



k,ita m +iV) k,ie> <!: +iv, 
h '<£>- k «B,„ 

j +■ j( V 

A ^ ^ B 

fcieo.+iv, fci'to„+!V, 




FIG. 6: Individual terms for the bare polarization in the or- 
dered phase. 



expression for the generalized polarization n m m+ / is de- 



picted in Fig. [6] and simplifies to 

TT — — \ " ■'2i r ,AA r BB 

ri m ,m+/ — , r / Jfel"fe ra "itra4 



^ km ^ km+l 



N 

k 

, i^BA/~<BA , (~<BB rtAA \ CX(\\ 

"T" Lj fem LT fcm+/ T u km u km+l ) ^ l JU 7 

where jk — — limi)^ 00 (i*/v A D) 2~2r*=i sm an< ^ solid 
lines denote the momentum-dependent lattice Green's 
functions G%> m [see Eqs. ([TUHT2]) ]. After substituting in 
the expressions for the Green's functions, the individual 
contributions to H m ,m+i become 

1 .^r AA r BB --(it, A- n B V B \ 

2_^,Jk (j km {j 'km+l ~ -^\ luJ m T A*d ^m) 



x(iuj m + in+ - S m+i ) 



F z (iu m+ i) _ F z (iuj m ) 
Z(iuJ m+ i) Z(iuj m ) 
Z 2 (iuj m ) - Z 2 (iuj m +i) ' 



\ " n 2 r BB r AA 

jy / t Jk Ksr km'~ r km+l 
k 



F z (iu m+ i) F z (iu> m ) 



Z(lUJ, 



m+l) 



Z{iu ri 



x(iuj m + iui + /if - ^ B l+l ) . 

Z Z (lUJ m ) - Z Z {lLU m+ l) 



AB 
km+l 



-2 /~iBA/~iBA 
Jk^km^km+l 



N ^ 

k 



(31) 



_1 Z(iuj m+ i)F z (iuj m+ i) - Z(iuj m )F z (iuj m ) 
2 Z 2 (iuj m ) - Z 2 (iuj m+ i) 

Hence, the full expression for Tl m ^ m+ i is 

F z (iuj m+ i) ^ F z (ioj m ) 
n _ 1 1 Z(iu m+ i) Z(iuj m ) 

m,m 2 i 2^{iuj m ) - Z 2 (iuj m+ i) 



x (iuj m + /if - E^)(iw m + ivi + fi d - T, m+l ) 

+ (iuj m + \x A - S^)(iw m + ivi + /if 

Z(iuj m+ i)F z (iuj m+ i) - Z(iu m )F z {iuj m ) 



yB \ 



Z 2 (ioj m ) - Z 2 (iu m+ i) | 

(32) 

Then, the expression for the current-current Green's 
function is obtained by substituting Eq. (f3"2"]l into Eq. (121)1) 
and analytically continuing the summation over Matsub- 
ara frequencies into contour integrations 



2m J 



x n(cD — i0 + , uj + ivi) — H(ui + i0 + , Q + iv{) 
+ n(o) — ivi, Cj — i0 + ) — H(Q — ii/i, uj + i0 + ) 



(33) 
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Here we have f(u>) = 1/[1 + exp(/3u))] is the fermi distri- 
bution function. The final step is to analytically continue 
from the bosonic Matsubara frequencies to the real axis 
(iv — > lu ± i0 + ). This produces our final result 



+00 



(2ni) 



dw [/(£)- /(« + «)] 



x Re{n(w - i0 + , u + oj + i0 + ) 
-n(w-iO + ,<D + w-iO + )}. (34) 



To make Eq. (|34| concrete, we substitute in the analytic 
continuation of Eq. (|3li|) to find the final expression for 
the optical conductivity (we set e 2 = 1): 



The final formalism we need to develop is for the dc 
transport properties. Starting from the expression for 
the optical conductivity in Eq. (|35|) we can calculate the 
dc conductivity by taking the zero frequency limit: 



a dc = lim cr(u!) 



(36) 



Air 2 



x Re<! 



F* z {Cj + uo) _ F z {Co) 
Z*(u! + ui) Z(uj) 
Z 2 {Cj)-[Z*{Cj + uj)] 2 



x (J<D + /Xd -Xr(<20P + w + /x£-£ A *(ui+u; 

+ \w + i4- [0 + w + /if - S s * (Cj + w)] 

Z*(^ + ^)F z *(^ + ^)-Z(^)f z (^) 

Z 2 (J>) - [Z*(J) + w)] 2 
F z (tD + o;) _ F z (d>) 
Z(Cj + w) ^(^) 



Z 2 (w)-Z 2 (tD + w) 
x ([u + nf -Y. b {Cj)][Cj + uj + -Y> a (Cj + uj)] 

+ [u + i4- S A (w)] [u> + uj + /tf - S S (w + w)] 
Z(w + w)F z (d) + w) - Z(C)F,(w) 1 



Z 2 {Cj) - Z 2 (Cu + oj) 



(35) 



The algebra is completely straightforward, but requires 
a careful use of l'Hopital's rule for determining some of 
the limits. After some lengthy algebra, we find that the 
final expression of the dc conductivity becomes 



G dc 



+00 

2 / dw 



4fM 



(37) 



with the exact many-body relaxation time t(u>) equal to 



"Re{[. 



u> ■ 



•i4 



■t4 



\ReZ(uj) ImZ(w) 



- 2Re[Z(u)F z (uj 
I 



gf>)]} [ ReF.M lmF z (uj) 
\ReZ(uj) lmZ(u) 



1] 



(38) 



For large frequencies the relaxation time approaches the 
asymptotic value 

1 2 

T °° = 4^ £/ 2 [nf (1 _ + n B (1 _ ; ( 39 ) 



this is a well-known anomaly on the infinite-dimensional 
hypercubic lattice^ due to the fact that the DOS never 
vanishes and at large frequencies the imaginary part 
of the self-energy is exponentially small, implying very 
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FIG. 7: (Color online) Exact many-body relaxation time at various T values for the CDW-ordered phase on a hypercubic 
lattice with (a) U — 0.5, (b) U = 1, (c) U = 1.5, and (d) U = 2.5. At high T in the Mott-insulator, the relaxation time goes 
to zero as u> 4 [this is hard to see in panel (d) because the quartic region occurs only for small frequencies and cannot be easily 
seen on this linear scale plot]. 



long lifetimes for the excitations. Note that the high- 
frequency limit of t(lo) actually diverges as T — > at 
half filling. This trend can be seen to develop in Fig. [3 
although we do not push the calculations too low in tem- 
perature due to accuracy issues with determining the sub- 
gap states. 

Starting from Eq. (|57|) we can also calculate the ther- 
mal transport. Since the system is at half-filling, the thcr- 
mopower vanishes clue to particlc-holc symmetry: the 
relaxation time in Eq. (|38[) is symmetric with respect to 
sublattice indices and is an even function of frequency at 
half- filling (Fig. [7]). The electronic contribution to the 
thermal conductivity K e is nonzero, and can be found in 
the standard fashion. It is expressed in terms of three 
different transport coefficients Ln, Lyi — L21 and L22 
as follows^ 



1 
T 



L22 — 



£12^21 
in 



In this notation, the dc conductivity satisfies 

Orfc = -^ll- 



(40) 



(41) 



The other transport coefficients can be calculated from 
the Jonson-Mahan theore m 20 ' 21 which says that there is 
a simple relation between these different transport coef- 
ficients, namely that they reproduce the so-called Mott- 
Thellung noninteracting form 2 ^, 



+00 



Lij — 



df(u) 



duj 



,i+i-2 



(42) 



where t(u>) is the exact many-body relaxation time de- 
fined in Eq. (38|) and plotted in Fig. 



III. NUMERICAL RESULTS 

We begin our discussion on transport properties in the 
ordered CDW phase by examining the optical conduc- 
tivity. In Fig. [5] we plot the temperature dependence of 
the optical conductivity for a dirty metal with U = 0.5. 
At high temperatures we see the expected behavior for a 
dirty metal — namely, there is a peak at low energy and 
a spread on the order of the metallic bandwidth. The 
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system does not have a low energy fermi liquid peak, 
because it is not a fermi liquid. Below the critical tem- 
perature for CDW order, the shape of the optical conduc- 
tivity changes significantly. Note how the spectral weight 
is shifted upward in frequency because the system is be- 
coming an insulator at low T. In particular, a sharp peak 
develops at u> = U which corresponds to the interband 
transitions from the lower band at uj < —U/2 to the up- 
per band at u> > U/2. We also see two additional peaks 
at lower frequencies. The higher of those peaks corre- 
sponds to transitions from the lower band to the subgap 
states above the chemical potential and from the subgap 
states below the chemical potential to the upper band 
and the lower one corresponds to the transitions between 
the subgap states below and above the chemical poten- 
tial. Both of these lower energy peaks must vanish as 
T — > because the subgap states disappear continuously 
at T = 0. Note that the frequency u> = U divides the 
spectra into two parts: to the right of this point the in- 
tensity of spectra increases as T decreases and to the left 
of this point the intensity decreases as T decreases which 
is similar to the isosbestic behavior of Mott insulators in 
the homogeneous phase, although we don't see the same 
kind of isosbestic behavior in the ordered-phase optical 
conductivity here. 

Results for U — 1 have a similar structure to those for 
U = 0.5, so we do not show them here. 

We next plot the optical conductivity for a near critical 
Mott insulator (U = 1.5) in Fig. O Here we see similar 
structure, with the peaks shifting to higher energy as T 
decreases, but the overall effect is not as large as in the 
metal, because this system would be an insulator even if 
there was no CDW order. Nevertheless, we still see the 
large peak develop with an edge at u> = U, and we see 
two low-energy peaks that have strong temperature de- 
pendence due to the types of transitions involving subgap 
states described above. 

Finally, we plot results for a moderate gap Mott insula- 
tor (U = 2.5) in Fig. [TTJ] The behavior here is essentially 
identical to what we saw at smaller values of U , except 
the effects are smaller, because the subgap states are very 
small for frequencies below where the Mott gap region ex- 
tends, so the overall effects are somewhat reduced. But 
all of the qualitative behavior remains. 

In order to complete our discussion of the dynamical 
response, we now describe the optical sum rule. In gen- 
eral, the sum rule for the optical conductivity is 



duia(uj) = —irK, 



(43) 
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FIG. 8: (Color online) Optical conductivity for U = 0.5 and 
various temperatures. Panel (a) is a linear scale and panel 
(b) is a logarithmic scale. 



netic energy is equal to 



m k 

= Ty][Z m F zm - 1] = T*y^ ^m^m = ^^^m^' 



= --J dto f{u) Ira \ a {uj)G a (uj), a = A,B, (44) 

and at T — 0, when S A = U and E s = 0, we immediately 
find 



K = -\J 



HI +f 2 
4 



(45) 



where K is the average kinetic energy (which is always 
nonpositive) . In the CDW-ordered phase the average ki- 



ln Fig. 1111 we plot the average kinetic energy both for 
the CDW and homogeneous solutions for different val- 
ues of U at T = 0. For small values of U (U < 0.648), 
we observe the anticipated behavior that the average ki- 
netic energy increases faster in the ordered phase than 
in the homogeneous phase. This is anticipated because 
the homogeneous phase has, on average, some neighbor- 
ing sites with no localized electrons, implying hopping 
is easier than in the ordered phase, where every hop in- 
volves a change in energy by U at T = because the 
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FIG. 9: (Color online) Optical conductivity for U = 1.5 and 
various temperatures. Panel (a) is a linear scale and panel 
(b) is a logarithmic scale. 



order parameter is uniform on each sublattice. Since it 
is more difficult to hop in the ordered phase, the kinetic 
energy increases relative to the homogeneous phase. For 
large values of U we find anomalous behavior, where the 
average kinetic energy is more negative in the ordered 
phase. There is no simple picture to explain how this 
occurs. In the homogeneous phase, as U increases, it 
becomes more difficult to hop because the doubly occu- 
pied states are being projected out of the system. This 
implies the average kinetic energy increases in the homo- 
geneous phase, but it does so faster than in the ordered 
phase. The subtle details of how the average kinetic en- 
ergy evolves with temperature are shown in Fig. 1121 The 
anomalous behavior for the temperature dependence of 
the average kinetic energy occurs for a finite range of T 
when U > 0.52. This is the "critical" U value where the 
DOS in the normal state changes its curvature from being 
negative at the chemical potential, as expected for a con- 
ventional metal, to positive in what is sometimes called 
an anomalous metal. In the region 0.52 < U < 0.648, the 
normal state DOS starts to develop a dip at the chemical 
potential, and for a finite temperature range, the anoma- 
lous behavior in the average kinetic energy occurs only 
for low temperatures. As U is increased further, we see 
the anomalous behavior occur for all T. These results 
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FIG. 10: (Color online) Optical conductivity for U = 2.5 and 
various temperatures. Panel (a) is a linear scale and panel 
(b) is a logarithmic scale. 



show that the spectral weight in the CDW phase shows 
a modest decrease for small U and a dramatic increase 
for large U at T = 0! This is somewhat unexpected, 
since the behavior is different than what is seen in say 
a BCS superconductor, where the gap formation reduces 
spectral weight at high frequencies, but the lost weight is 
restored in a zero frequency Drude peak. For the CDW 
ordered phase, no zero frequency delta function appears. 
The spectral weight loss is small for small U, but the gain 
can become significant for large U. 

Next we examine the dc transport. The temperature 
dependence of the dc and thermal conductivity are plot- 
ted in Figs. [T51 and [T4l respectively, where we plot both 
the CDW solution and the homogeneous solution extrap- 
olated into the CDW region. At low temperatures, due 
to the factor —df{uj)/dLj^ the main contributions to the 
dc transport come from the narrow region of width 4T 
around the chemical potential (the so-called fermi win- 
dow) . For the Falicov-Kimball model at half filling in the 
homogeneous phase (T > T c ) the DOS, Green's functions 
and self-energies do not depend on temperature and, as 
a result, the temperature dependence of the dc transport 
is determined solely by the shape of the relaxation time 
in Eq. (|3S[) close to the chemical potential. For small U 
values the relaxation time t(ui) is flat [Fig. [7] (a)] and, as 
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a result, the dc conductivity for the homogeneous phase 
is essentially a constant for low T. At U = V2 the Mott 
insulator forms. For larger U values, one might expect 
to see exponentially activated transport, but that does 
not occur on the hypercubic lattice, because the system 
only possesses a pseudogap. Even though the DOS expo- 
nentially decreases in the gap region, the lifetime of the 
excitations is exponentially long, and t(oj) behaves like 
uj a for low energies^. This produces a quartic depen- 
dence of the dc conductivity on T, and a higher power 
law for the thermal conductivity. 

In the CDW phase (T < T c ), the CDW gap is filled 
by subgap states at finite T, which lead to a less severe 
modification of the exponentially activated transport at 
low T. But it is only the subgap states within the fcrmi 
window that affect the transport, so the modification is 
not quite as severe as one might have naively guessed. 
Note the small wiggles in the solid lines at low T. These 
occur due to the evolution of the subgap states. The T 
dependence of the dc transport always shows a marked 
kink at T c with the conductivities sharply suppressed as 
the CDW gap forms. In the Mott insulator, the trans- 
port changes from power law in T to exponential acti- 
vation (suitably modified by the subgap states). The 
thermal conductivity displays similar features, as shown 
in Fig. [H 



IV. CONCLUSIONS 

In this work we have developed the formalism to calcu- 
late transport properties of CDW-ordered phases within 
DMFT. Since the dc charge and heat transport and the 
optical conductivity continue to have no vertex correc- 
tions, even in the ordered phase, the calculations reduce 
to a careful evaluation of the bare Feynman diagrams 
with a sublattice index introduced by the order. 

As the system orders into a CDW state, the DOS de- 
velops a gap with a sharp singularity in the DOS at the 



0.00 0.01 0.02 0.03 
Temperature [t*] 



0.04 




-0.247 

>> 

is 

g -0.248 

CD 
C 

-0.249 

0.00 0.01 0.02 0.03 0.04 0.05 0.06 
Temperature [t*] 

-0.210 

>, 

iu -0.215 



| -0.220 

-0.225 
0.00 




0.02 0.04 
Temperature [t*] 



FIG. 12: Temperature dependence of the average kinetic en- 
ergy for different values of U: (a) U = 0.5; (b) U = 0.52; 
(c) U = 0.645; (d) U = 1. The solid line corresponds to the 
CDW phase and the dotted line corresponds to the homoge- 
neous solution. 



band edge when T = 0. The gap at T = is always equal 
to U. As the temperature increases, but still below T c , 
we see a significant development and evolution of sub- 
gap states within the gap region. This gap region where 
subgap states develop, appears to lie within the extent 
of the normal-state DOS — in other words, in the Mott 
insulator, we do not see subgap states develop within the 
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FIG. 13: dc conductivity for (a) U = 0.5 (T c « 0.034), (b) U = 1 (T c « 0.0615), (c) t/ = 1.5 (T c « 0.075) and (d) (7 = 2.5 
(T c w 0.072) as a function of temperature. The solid line denotes the CDW ordered phase and the dashed line denotes the 
homogeneous one. Results are presented on a linear scale (left) and logarithm of dc conductivity vs inverse temperature (right). 



region that corresponds to the Mott gap in the normal 
state. We verify the accuracy of the DOS calculations by 
calculating the zeroth, first, and second moment of the 
local DOS on each sublattice and we find they agree with 
exact results to essentially machine accuracy. 



The optical conductivity has a significant rearrange- 
ment of states within the ordered phase, which can be 
understood by examining the different kinds of processes 
that take place within an optical transition — namely that 
we move from an occupied to an unoccupied state. Be- 
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FIG. 14: Electronic contribution to the thermal conductivity for (a) U = 0.5 (T c w 0.034), and (b) U = 2.5 (T c w 0.072) as 
a function of temperature. The solid line denotes the CDW ordered phase and the dashed line denotes the homogeneous one. 
Results are presented on a linear scale (left) and logarithm of thermal conductivity vs inverse temperature (right). 



cause there are many different bands that are present at 
finite T, this leads to significant structure in the optical 
conductivity. In particular, the singularity in the DOS 
leads to a large asymmetric peak centered around U in 
the response function. The total spectral weight is gov- 
erned by the average kinetic energy due to the optical 
sum rule. While a naive expectation would say the aver- 
age kinetic energy increases when the ordering is turned 
on (i. e., it becomes less negative with a smaller magni- 
tude) because the ordering blocks hopping between the 
sublattices, we find that is true only for small U. For 
small U the kinetic energy shows a modest increase, so 
some spectral weight is lost due to the ordering. For 
larger U the kinetic energy shows a significant reduction 
(i. e., the magnitude increases as the average kinetic 
energy becomes more negative) so the spectral weight 
increases when the ordered phase is entered, and that in- 
crease can become quite substantial as U becomes large. 

Finally, we also examined the dc transport. Since we 
are at half filling, one can show the thermopower vanishes 
due to particle-hole symmetry even in the presence of 
CDW order. Hence we can only examine the charge and 
heat conductivities. We find that the CDW order sup- 
presses both of these, but because of the subgap states 
and their complicated evolution with temperature, the dc 
response does not obey any simple functional form at low 



T. Instead, we often see significant wiggles in the conduc- 
tivities. In the Mott-insulating phase, the conductivity 
should go from a power-law-like behavior to exponential 
activation. We see such a trend start to develop, but 
we cannot accurately quantify this because we cannot go 
down far enough in temperature in the CDW phase be- 
fore we run into issues with accuracy of the calculations. 

This work shows that there is rich and interesting be- 
havior that occurs in the transport as CDW order sets 
in. In future work, we will examine Raman scattering, 
where vertex correction effects arc present and inelastic 
X-ray scattering, where interesting phenomena is likely 
to occur when the photon transfers momentum equal to 
the ordering wavevector. 
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